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1. Introduction 

In this paper we review briefly the generalized quantum dynamics 1 ' 2 constructed on 
a phase space of local noncommuting fields. We show that the equilibrium conditions on 
the microcanonical entropy imply that the system decomposes thermodynamically to a 
sequence of adiabatically independent subsystems, each with its own temperature. There 
is an equipartition theorem for the phase space variables of the system generated by the 
linear combination of conserved quantities associated with each of these independent ther- 
modynamic modes. 

We start with a review of our basic framework. Generalized quantum dynamics 1 ' 2 is 
an analytic mechanics on a symplectic set of operator valued variables, forming an operator 
valued phase space S. These variables are defined as the set of linear transformations^ on 
an underlying real, complex, or quaternionic Hilbert space (Hilbert module), for which the 
postulates of a real, complex, or quaternionic quantum mechanics are satisfied 2 " 6 . The 
dynamical (generalized Heisenberg) evolution, or flow, of this phase space is generated by 
the total trace Hamiltonian H = TrH, where for any operator O we have 

O = TrO = ReTr(-l) F 

= ReJ2(n\(-l) F 0\n), ( L1 ) 

n 

H is a function of the operators {q r (t)}, {p r (t)}, r = 1,2, ...,N (realized as a sum of 
monomials, or a limit of a sequence of such sums; in the general case of local noncommuting 
fields, the index r contains continuous variables), and (— 1) F is a grading operator with 
eigenvalue 1(— 1) for states in the boson (fermion) sector of the Hilbert space. Operators 
are called bosonic or fermionic in type if they commute or anticommute, respectively, with 
(— 1) F ; for each r, p r and q r are of the same type. 

The variation of a total trace functional with respect to some operator is defined with 
the help of the cyclic property of the Tr operation. The variation of any monomial O 
consists of terms of the form Ol8x t Or, for x r one of the {q r }, {Pr}, which, under the Tr 
operation, can be brought to the form 



50 = 5TrO = ±TrO R L Sx r , 



so that sums and limits of sums of such monomials permit the construction of 

50 = Tr^^-<kv, (1.2) 

r 

uniquely defining 50/5x r . 

Assuming the existence of a total trace Lagrangian 1 ' 2 L = ~L({q r }, {q r }), the variation 
of the total trace action 

/oo 
L({q r },{q r })dt (1.3) 
-oo 



t In general, local (noncommuting) quantum fields. 
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results in the operator Euler-Lagrange equations 

oL d 5L 



0. (1.4) 



8q r dt 5q r 

As in classical mechanics, the total trace Hamiltonian is defined as a Legendre transform, 

H = Tr ^ p r q r — L, (1.5) 

r 

where 

Pr=J~- 1-6 



It then follows from (1.4) that 



T— = -Pr T— = Cr?r, (1-7) 



where e r = 1(— 1) according to whether p r ,q r are of bosonic (fermionic) type. 
Defining the generalized Poisson bracket 

{AB}-TrVe (——-— — } (18a) 
^ r \Sq r Sp r Sq r Sp r ) ' 

one sees that 

dA dA . ,„ 

1F = 1F + {A ' H} - (L86) 

Conversely, if we define 

x s {r)) = Tr{r)x s ), (1.9a) 

for rj an arbitrary, constant operator (of the same type as x s , which denotes here q s or p s ), 
then 

dx ^) - rj^x^^f tes (rj) oH oH ox s (77) \ 

r \ 5(7r 0£> r 5(7^ Sp r J 

and comparing the coefficients of 77 on both sides, one obtains the Hamilton equations (1.7) 
as a consequence of the Poisson bracket relation (1.86). 

The Jacobi identity is satisfied by the Poisson bracket (1.8a), 7 and hence the total 
trace functionals have many of the properties of the corresponding quantities in classical 
mechanics. 8 In particular, canonical transformations take the form 

5x s (r/) = {x s (r/),G}, (1.10a) 

which implies that 

5p r = -—, 5q r = e r — 7 (1.106) 
oq r op r 



with the generator G any total trace functional constructed from the operator phase space 
variables. Time evolution then corresponds to the special case G = Hdt. 

It has recently been shown by Adler and Millard 9 that a canonical ensemble can be 
constructed on the phase space S, reflecting the equilibrium properties of a system of many 
degrees of freedom. Since the operator 

C = ^2(e r q r p r -p r q r ) 

(1-11) 

r,B r,F 

where the sums are over bosonic and fermionic pairs, respectively, is conserved under the 
evolution (1.7) induced by the total trace Hamiltonian, the canonical ensemble must be 
constructed taking this constraint into account. This is done by including in the canonical 
exponent the conserved quantity TrA(7, for some given constant anti-hermitian operator 
A. 

In the general case, in the presence of the fermionic sector, the graded trace of the 
Hamiltonian is not bounded from below, and the partition function may be divergent. 
When the equations of motion induced by the Lagrangian L coincide with those induced 
by the ungraded total trace of the same Lagrangian, 

L = ReTrL, 

without the factor (— 1) F , the corresponding ungraded total trace Hamiltonian H is con- 
served; it may therefore be included as a constraint functional in the canonical ensemble, 

along with the new conserved quantity TrXC (see Appendices and C of ref. 9), where 



C = ^2[q r ,Pr] 

(1.12) 

r,B r,F 

It was argued that the Ward identities derived from the canonical ensemble imply that 

A and A are functionally related, so that they may be diagonalized in the same basis 
(Appendix F of ref. 9). It was then shown that, since the ensemble averages depend only 
on A and (— 1) F , the ensemble average of any operator must commute with these operators. 
Since the ensemble averaged operator (C)av is anti-self-adjoint, if one furthermore assumes 
it is completely degenerate (with eigenvalue i e ffh), the ensemble average of the theory then 
reduces to the usual complex quantum field theory. 

As discussed in detail in ref. 10, the phase space volume associated with the micro- 
canonical ensemble can be written as 

T(E, E, v, v) = I ' dfi6(E -H)S(E - H) 

Sr " (1-13) 

II t(»L-(n\(-l) F C\m) A )6(^ m -(n\C\m) A ), 

n<m,A 
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where we have taken into account the possible algebraic structure of the matrix elements 
of the operators with the index A, which takes the values 0, 1 for the complex Hilbert 
space, 0, 1, 2, 3 for the quaternionic Hilbert space, and just the value for real Hilbert 
space. The invariant phase space measure is defined by 

r,m,n 

where redundant factors are omitted according to adjointness conditions. We have, fur- 
thermore, used the abbreviations v = {v^m} an d v = {^nm}- The entropy associated with 
this ensemble is given by 

S mic (E, E, i>, v) = In F(E, E, i>, v). (1.15) 

It was argued in ref. 10 that a large system can be decomposed into a part within a 
certain (large) region of the measure space, which we denote as 6, corresponding to what 
we shall consider as a bath, in the sense of statistical mechanics, and another (small) part 
which we shall denote as s, corresponding to what we shall consider as a subsystem. It was 
then argued 10 that the phase space volume can be well approximated by 

T{E, E, i>,v) = J dE s dE s (dv s )(dv s ) T b (E -E S ,E- E s , v-v 8 ^v- P s ) T S (E S , E s , z> s , i> s ). 

(1.16) 

Defining the variables 

S = {Si} = {E,E,v,i>}, (1.17) 

it was shown 10 that the equilibrium conditions which follow from the assumption that 
there is a maximum in the integrand of Eq. (1.16) (which dominates the integral in the 
limit of a large number of degrees of freedom) result in the set of equalities 

1 ar, -K)lf=P^-«££(s-0lf. a-") 



where S corresponds to the total quantities belonging to the full system. It was then shown 
that the canonical ensemble obtained by Adler and Millard 1 ' 2 , 

p = Z- 1 exp-{rH + fH + TVAC + frAC'}, (1.19) 

where 

Z = J d/JL exp -{rH + fH + TrXC + frxd}, (1.20) 

follows in a straightforward way. The quantities r, f and the matrices (real, complex, or 
quaternionic) A, A are the equilibrium parameters defined by the values of the members 
of (1.18) for each of the £'s 10 ; they therefore correspond to temperatures precisely as they 



5 



emerge in conventional statistical mechanics. We remark that Ingarden 11 has studied a 
similar generalization of temperature in the framework of the statistical mechanics associ- 
ated with problems of optical pumping (in the diagonal form which we shall discuss in the 
next section). 

Replacing the operators and trace functionals in (1.20) by integrals over ^-functions, 
the partition function can be rewritten as 10 

Z = J ' dEdEidu^d^e^^'^exp-irE + tE + Tr\i> + TrU}. (1.21) 

By studying the dispersions of the variables in the canonical ensemble, it was found 10 
that the second derivative matrix of the micro canonical entropy is negative definite, i.e., 
that 

< L22 > 

In the following we use the fact that this matrix is real symmetric to diagonalize it, 
and in this way to construct a set of dynamical generators over which the total entropy 
decomposes in a neighborhood Co of the maximum entropy point. 

2. Diagonal form of the second variation of the entropy 

The negative definite matrix Eq. (1.22) 

p - mtc (2 1) 

13 " diM 3 ( 1 

is symmetric and can therefore be diagonalized by an orthogonal transformation. Let 
(orthogonal) be such that, in the neighborhood Co, 

y]flfc«fl£jAj = faidjjt), (2.2) 

ij 

where the elements g^(£) on the right hand side are the negative eigenvalues. Now, let us 
define, using these constant coefficients, 

efc = y]flfcj&, (2.3a) 



and hence 



It then follows that, in Co, 



& = y^Qfcjefc- (2.36) 

k 



d 2 S ^ d 2 S 



(2.4) 



1-3 
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Since the crossed derivatives of S vanish, S must be a sum of functions that depend on 
each of the {e k } separately, i.e., 

S = J2Sk(e k ). (2.5) 

k 

The entropy is therefore additive (in Co) over diagonal "thermodynamic modes". 
The equilibrium parameters defined in the previous section, 

OS 

Xj = W = {r ' f ' A ' A} ' (2 ' 6) 

may be transformed in the same way, i.e., 



\- d dS 

dS k (e k ) __ 1 



8e k T k " Tl de k ' T fc 2 C k 
are positive, and by (2.7) — (2.9) are given by 



(2.7) 



de k T fc ' 

giving the diagonal temperatures (of the type considered by Ingarden 11 ). 
We remark that, according to (2.2) and (2.4), 

d 2 S 

84= dt< °' {2 ' 

so that the "specific heats" , entering as 

d 1 1 dT k 11 



(2.9) 



C k = -^d k . (2.10) 



rp2 



3. Equipartition 

We now consider linear combinations of the dynamical quantities 

Hi = {H, H, C, C} 

of the same form as the linear combinations of the parameters {^} which are their equi- 
librium values, 

£fc = ^2 a kiHi, (3.1) 



the effective "energies" associated with the thermodynamic modes. Since the determinant 
of the matrix a is unity, the microcanonical phase space integral (1.13) can be written as 

r(£) = ! dn\[8{e k -e k ); (3.2) 
J k 

since, however, as we have shown in Section 2, 

lnr(£) = S( ei ,e 2 ,...) 

k 

it follows that the phase space volume factorizes on the diagonal parameters 



(3.3) 



k 

= r( ei ,e 2 ,...) = n r *( e *)- 



(3.4) 



One can show that the free energy also becomes additive 12 . 
Let us now consider the microcanonical average 

(Vr-p-) = TT t " r / dnT[5(e e - e e )q r -^, (3.5) 

5q s ' r(ei,e 2 , ...) J Y oq s 

where {q r } are the canonical coordinates (fields) of the phase space. We now write the 
right hand side of (3.5) identically as 

w ~ vIT^r - / d ^ QrJ-(£k -e fc ), 

r(ei, e 2 , . . .) Y dei J {Sj<ej} 5q s 

replacing the 5-functions by derivatives of the parameters of boundary step functions; 
adding the constant e k does not affect the result. Integrating by parts in the integration 
over phase space, we obtain 

(Qr^) = Tr t f {Y~Qr(£k - e fc ) - S rs (e k - e fc )}. (3.6) 

x 5q s ' r(ei,e 2 ,...) Y de e J e . <e . 1 Sq s 

The first term vanishes on the boundary, and we therefore have 

<*£> = -f7r^— (3-7) 



5q s / r( ei ,e 2 ,...) Y^, J {sj<ej} 

The derivative with respect to e k in the product of derivatives vanishes when it differ- 
entiates the upper bound; its contribution is only from the integrand, resulting in a factor 
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—1. The other derivatives act only on the upper limits. The product then results in the 
restricted measure 

/ dfji(ee = ei V k) = E fc , (3.8) 

J e k <e k 

which can be rewritten as 

S fc = / dnj[8{e t -e t ). (3.9) 

According to (3.2), 

^ = r( ei ,e 2 ,...). (3.10) 

(irp) = =r, ^ ~^ k . (3.11) 
5q 8 ' r(ei,e 2 ,...) 

We now use the factorization of T(ei,e2, . . .) in Co to derive a relation between and 
the additive entropies. In the limit of a large number of degrees of freedom, the leading 
edge of the integral defining £& dominates the integral 13 , so we may formally extrapolate, 
as a model, the quadratic form (and associated factorization) valid in Co- From (3.4), and 
(3.10) it then follows that 

ay 



We therefore have 



One may integrate this equation to obtain 

S fc = I" T k (e' k )de' k \[T t {et) +G(e t ,l^ k). (3.13) 

l±k 

Since e k cannot be — oo (the functional H is contained linearly and its positive values are 
assumed to dominate for large values of the phase space variables), the first term on the 
right hand side of (3.13), along with T, k must vanish as e k — > — oo, and hence G must be 
zero. 

We therefore obtain 

/efc 
de'k eSfe(6fc) II e ' ( 3 - 14 ) 
l±k 

so that 



J ek de' k e s ^ 



r (ei ,e 2 ,...) - e*c) (3J5) 
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With the leading approximation for a large number of degrees of freedom 13 , 

In / de' k e Sfc(e ' fe) ~ lne Sfe(efe) , 



we conclude that 



((ZrT ) = —8 rs Tk. 



(3.16) 



We finally make some remarks on the flows generated by which, for clarity, we 
recast to the form (summed on nm) 



The Poisson bracket (1.8a) then contains a term which is the t-derivative, by (1.96), 
but there are additional terms of general type (1.106). In Ref. 10, it is shown that the 

terms in (3.17) which contain C, C induce transformations on phase space which are 
commutators with a^, in the boson sector, and with in the fermion sector, but anti- 
commutator with 5/s in the fermionic sector. Hence the elements of the diagonalization 
transformation act as connection forms under evolution generated by the effective mode 
energy functionals. Further discussion and application of these results will be given in Ref. 
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